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We discuss various sources of azimuthal correlations in relativistic heavy-ion collisions. The inte- 
gral measure <3> is applied to quantify the correlations. We first consider separately the correlations 
caused by the elliptic flow, resonance decays, jets and transverse momentum conservation. An ef- 
fect of randomly lost particles is also discussed. Using the PYTHIA and HIJING event generators 
we produce a sample of events which mimic experimental data. By means of kinematic cuts and 
particle's selection criteria, the data are analyzed to identify a dominant source of correlations. 
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I. INTRODUCTION 

Particles produced in relativistic heavy-ion collisions 
are correlated in azimuthal angle due to various mecha- 
nisms. One mentions here extensively studied jets and 
minijets resulting from (semi-)hard parton-parton scat- 
tering and collective flow due to the cylindrically asym- 
metric pressure gradients, see the review articles [1] and 
[2], respectively. More exotic sources of correlations are 
also possible. As argued in [3], the plasma instabilities, 
which occur at an early stage of collisions, can generate 
the azimuthal fluctuations. Except the dynamically in- 
teresting mechanisms, there are also rather trivial effects 
caused by decays of hadronic resonances or by energy- 
momentum conservation. 

There is a variety of methods designed to study fluc- 
tuations on event-by-event basis. In particular, the so- 
called measure $ proposed in [3] was used to measure 
the transverse momentum [5J [5] and electric charge fluc- 
tuations [TJ. The measure proved to be very sensitive to 
dynamical correlations and it was suggested to apply it 
to study azimuthal ones [5]. Such an analysis is under- 
way using experimental data accumulated by the NA49 
and NA61 Collaborations and some preliminary results 
are already published [9]. The aim of this paper is to 
present model simulations to be used in interpretation of 
the experimental data. The fact that the measure $ is 
sensitive to correlations of various origin is advantage and 
disadvantage at the same time, as it is difficult to disen- 
tangle different contributions. Therefore, we model the 
azimuthal correlations driven by several processes and 
we look how the correlations show up when quantified by 
the measure We first consider separately in terms of 
toy models, the elliptic flow, resonance decays, jets and 
transverse momentum conservation. An effect of ran- 
domly lost particles is also examined. Then, we analyze 
the data provided by the PYTHIA and HIJING event 
generators showing how to identify the main sources of 



correlations by applying kinematic cuts and particle's se- 
lection criteria. 

A magnitude of $ of azimuthal correlations measured 
in Pb-Pb collisions at 158A GeV is of order of 0.01 ra- 
dian [9] . Since the effects, which are theoretically studied 
here, often exceed this order, a realistic model should in- 
clude all the effects, and consequently, the model has to 
be rather complex. For this reason we do not attempt to 
compare our model calculations to the preliminary data 
[5] . There is also an additional reason. The measurement 
[5] has been performed in an acceptance window which 
is not only limited in particle's rapidity and transverse 
momentum but it is also nonuniform in azimuthal an- 
gle. The acceptance, which is characteristic for the NA49 
detector, has to be properly included in model calcula- 
tions for a quantitative comparison with the experimental 
data. It is beyond the scope of this study. Our aim here is 
to understand how the collective flow, resonance decays, 
jets and transverse momentum conservation contribute 
to azimuthal correlations and how the contributions can 
be disentangled. 



II. MEASURE $ 

Let us first introduce the correlation measure <I>. One 

def 

defines the variable z = x — x, where a; is a single par- 
ticle's characteristics such as the particle transverse mo- 
mentum, electric charge or azimuthal angle. The over- 
line denotes averaging over a single particle inclusive dis- 
tribution. In the subsequent sections, x will be identi- 
fied with the particle azimuthal angle (f> and the fluctu- 
ation measure will be denoted as The event vari- 
able Z, which is a multiparticle analog of z, is defined 

as Z = J2iLi( x i ~ ^0) where the summation runs over 
particles from a given event. By construction, (Z) = 0, 
where (...) represents averaging over events (collisions). 
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The measure $ is finally defined as 



def 



(AO 



(1) 



It is evident that $ = 0, when no inter-particle corre- 
lations are present. The measure also possesses a less 
trivial property - it is independent of the distribution of 
the number of particle sources if the sources are identi- 
cal and independent from each other. Thus, the measure 
$ is 'blind' to the impact parameter variation as long as 
the 'physics' does not change with the collision centrality. 
In particular, 3> is independent of the impact parameter 
if the nucleus-nucleus collision is a simple superposition 
of nucleon-nuclcon interactions. In the following sections 
we discuss how various mechanisms responsible for az- 
imuthal correlations contribute to Then, using the 
event generators we show how the dominant contribu- 
tions can be identified. 



III. COLLECTIVE FLOW 

Particles produced in relativistic heavy-ion collisions 
reveal a collective behavior which is naturally described 
in terms of hydrodynamics where the collective flow is 
caused by pressure gradients. The collective flow as quan- 
tified by the measure $0 was studied in [5] . Here we only 
recapitulate the results derived in [8]. 

Since the inclusive azimuthal distribution is flat, <h = tt 



and 



A2 - 42 



for <fi £ [0, 27r], and thus z 2 = \in 2 . A single 



particle azimuthal distribution in a given event is 



P{cj>) 



I 

2/7 



2 v n cos(n((p - <j>n)) 



(2) 



where < <f> < 27r; (f>ji is the reaction plane angle and 
v n denotes an amplitude of the n— th Fourier harmonics. 
The N— particle distribution is assumed to be a product 
on A" distributions ([2]) multiplied by a multiplicity dis- 
tribution. Consequently, the collective flow is the only 
source of azimuthal correlations in the system. Averag- 
ing over particles is performed by integrating over <pi with 
i = 1 , 2, . . . A" and averaging over events is achieved by 
integrating over (j)R and summing over N. The distribu- 
tion of reaction plane angle is obviously flat. Thus, one 
finds the measure $ of azimuthal correlations caused by 
the flow as 



(3) 



where (N m ) is the m— th moment of multiplicity distri- 
bution and 



We have first verified the effect of second Fourier co- 
efficient V2 on $0 by Monte Carlo simulations. For this 
purpose we have generated events of particle multiplic- 
ity given by either Poisson or Negative Binomial (NB) 
distribution. The latter is defined as 



r(N + k) (N) N k k 

T(N + l)T(k) ((jv)+fc) 



Pn = 



N+k • 



(5) 



where T(k) is the Gamma function, which for positive 
integer arguments equals T(k) = (k— I)!; the parameter k 
can be expressed through the variance of the distribution 
Var(AT) = (N 2 ) - (N) 2 and the average value (AT) as 



k = 



(N)'< 



Var(A0 - (AT) ' 



(6) 



The parameter k is cho sen in such a way in our all sim- 
ulations that y / Var(AT) = (N)/2. Then, the multiplicity 
distribution approximately obeys the Wroblewski's for- 
mula [10 which is known to hold for proton-proton in- 
teractions in a wide collision energy range. The simula- 
tions are performed for both the Negative Binomial and 
Poisson distributions as the former distribution is much 
broader than the latter one for (N) ^> 1. We note here 
that the width of multiplicity distributions in relativis- 
tic heavy-ion collisions strongly depends on centrality se- 
lection criteria. Thus, it is important to see how the 
correlation signal changes with the width of multiplicity 
distribution. 

The azimuthal angle of each particle has been gener- 
ated from the distribution 



P{4>) 



i 

2n 



l + 2i; 2 cos(2( 



0) , 



(7) 



where < (j) < 27r; the reaction plane angle ipR of a given 
event has been generated from the flat distribution. The 
results of our simulations are shown in Fig. [T] for both 
the Poisson (left panel) and NB (right panel) multiplicity 
distributions. As seen, the analytical formula (J3j) works 
perfectly well. 

There are large (~ 40%) event-by-event fluctuations of 
i>2 observed [2] at BNL RHIC. The v 2 fluctuations are 
dominated by the fluctuations of eccentricity of the over- 
lap region of colliding nuclei, see e.g. [TT]. We have in- 
troduced the V2 fluctuations in our simulations in the fol- 
lowing way. For each event the value of v 2 has been gen- 
erated from the Gaussian distribution of the dispersion 
a V2 . The fluctuations have been restricted to vary within 
2a V2 around the mean (1*2) that is (^2) — 2a V2 < v 2 < 
(v 2 ) + 2cr.„ 2 . Then, v 2 remains positive unless <r V2 /(v 2 ) 
exceeds 0.5. 

In Fig. [2] we demonstrate the effect of flow fluctua- 
tions relative the effect of flow. Specifically, we show the 
difference of the correlation measures computed for the 
fluctuating v 2 and fixed v 2 = (1*2) • The particle multi- 
plicity has been generated from the NB distribution with 
(A 7 ) = 400. As seen, the flow fluctuations of relative 
magnitude of ~ 40% noticeably increase the value of $0 
if (v 2 ) is not too small. 
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FIG. 1: (Color online) as a function of the second Fourier coefficient V2 for the Poisson (left panel) and NB (right panel) 
multiplicity distributions. The lines represent the analytical formula j3J. The insets show for small values of vi. 



IV. RESONANCE DECAYS 

Let us start the discussion of effects of resonance decays 
with the toy model where all produced particles come 
from heavy resonances which have vanishing transverse 
velocity and decay back to back into pairs of particles. 
The particle multiplicity is arbitrary but fixed even num- 
ber. Then, as shown in Appendix, we have 



1 - y/2 

$^ = j^ix « -0.531. 

V6 



(8) 



When only a fraction / of all produced particles comes 
from the back-to-back decays of resonances while the re- 
maining particles are produced independently from each 
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FIG. 2: (Color online) The difference of $0 computed for the 
fluctuating and fixed v 2 . 



other, the calculations presented in Appendix lead to 



$a = 



V6 7 



(9) 



As seen, for / = the formula ^ gives, as expected, 
$ = and for / = 1 we get the value 

We have checked the formula Q by Monte Carlo simu- 
lations and then we have considered the model where the 
particle multiplicity is not fixed but it is given by either 
Poisson or NB distribution with the average multiplicity 
equal 50, 400, or 700 particles. For a given fraction / 
of particles coming from the back-to-back decays of res- 
onances, the number of particles coming from the decays 
in the event of multiplicity N has been the even number 
which is the nearest to fN. The measure <E>^ as a func- 
tion of / is shown in Fig. [3] As seen, the formula ^ still 
works very well. 

When a resonance, which is at rest, decays back to 
back, the difference of azimuthal angles of the decay 
products is Ac/> — tt. When the resonance has a finite 
velocity, the difference of azimuthal angles of the decay 
products is smaller than tt. When the resonance's kinetic 
energy is much larger than the energy released in its de- 
cay, Ac/) is zero. Therefore, we consider a model where 
a fraction of particles comes from the resonance decays 
and the particles are emitted in pairs with the difference 
of their azimuthal angles Ac/) varying from to tt. We 
first assume that all particles are emitted in pairs and 
the relative azimuthal angle of two correlated particles 
equals Ac/>. As shown in Appendix, we then have 



9,, = 



A«/>7r + -(A</>)2--^. 



(10) 



As seen that changes its sign from positive to negative 
with growing Ac/>; $0 vanishes when 

7^ 



Ac/> = 7T 1 



1.328 



(11) 
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FIG. 3: (Color online) $0 as a function of fraction of particles coming from the back-to-back resonance decays. The particle 
multiplicity is distributed according to the Poisson (left panel) or NB (right panel) distribution. The dashed and solid lines 
represent the analytical formulas Q and qQp , respectively. 



and for = 7r we deal with the model described by 
the formula ([8| . Further on, we have considered a model 
where only a fraction / of particles is emitted in corre- 
lated pairs. Then, as explained in Appendix, Eq. ( 10 ) 
gets the form 



V3' 



(12) 



As previously, $0 changes its sign and $0 = for A<j) 
given by Eq. |TT] ) . 

In Fig. [4] we compare the formula ( 12 1 with the results 
of Monte Carlo simulation of <J>0 as a function of A<j>. The 
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FIG. 4: (Color online) $0 as a function of the correlation 
angle A<f> for varying fraction / of particles emitted in pairs. 
The particle multiplicity is given by the NB distribution. 



fraction of particles emitted in pairs equals 0.3, 0.5, 0.7 
or 1.0. The remaining particles, which are not emitted 
in pairs, carry no correlations. The particle multiplicity 
is generated from the NB distribution with ^/Var(TV) = 
(N)/2 — 50. As seen, the formula (12) works perfectly 
well. 

We next discuss the combined effect of resonance de- 
cays and elliptic flow. The multiplicity of events is gener- 
ated from the Poisson or NB binomial distribution. For 
each event 30% of particles is assumed to originate from 
heavy resonances which decay back to back into pairs of 
particles. Neither resonances nor their decay products 
experience any flow, but the remaining 70% of particles 
manifest the collective elliptic flow according to Eq. ([7]). 
The results of the simulation are shown in Fig. [5] When 
the particle multiplicity or V2 are sufficiently small, the 
effects of resonance decays dominates and $^ is nega- 
tive. It becomes positive when the effect elliptic flow 
takes over. 



V. 'DUETS' 

We call a dijet the two groups (jets) of particles flying 
in exactly opposite directions. Particles from each jet are 
distributed within a cone of the azimutal angle a<f). We 
have considered the dijets of 2+2, 5+5 and 10+10 parti- 
cles. The total particle multiplicity is correspondingly 4, 
10, 20 as there is exactly one dijet per event and there 
are no other particles. The results of dijet simulation are 
shown in Fig. [6] 

We have here two sources of azimuthal correlations 
which counteract each other. As we already know, the 
back-to-back emission of particles generates negative cor- 
relations while the collinear emission leads to positive 
ones. When the particle's multiplicity of dijets is suf- 
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FIG. 5: (Color online) $0 resulting from the combined effect of elliptic flow and resonance decays. The particle multiplicity is 
distributed according to the Poisson (left panel) or NB (right panel) distribution. 



ficiently high and a<p is sufficiently small, the effect of 
collincar emission wins and is positive. 



VI. MOMENTUM CONSERVATION 

The momentum conservation obviously leads to inter- 
particle correlations. We have studied the effect on <J>^, 
generating the sets of particles of multiplicity N. The az- 
imuthal angle distribution of a single particle is assumed 
to be fiat while the transverse momentum distribution is 
chosen in the form 



P{ PT ) = p 2 p T e 



-Ppi 



(13) 



with the slope parameter (3 1 = 200 MeV. For each par- 
ticle the x and y components of its momentum have been 
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FIG. 6: (Color online) $0 as a function of the jet opening 
angle a<j> f° r several numbers of particles in a jet. 



computed as p x = pxcoscj) and p y — p^sm/p. To make the 
total transverse momentum of N particles vanish, the x 
and y component of momentum of each particle has been 
shifted as 



1 

Px -> Px - -T7 



N 

v 51 

1=1 



Px , Py^Py 



1 N 

»=1 



(14) 



The simulation showing the effect of transverse mo- 
mentum conservation is illustrated in Fig. [TJ The parti- 
cle multiplicity has been generated according to NB dis- 
tribution with y/Vai(N) — (N)/2. As seen, the effect 
of momentum conservation is sizable and it survives to 
large multiplicities. 

In real experiments only a fraction of all produced par- 
ticles is observed due to a finite detector efficiency and 
acceptance. We model the effect of detector efficiency 
by randomly loosing particles independently of their az- 
imuthal angle. In Fig. [8] we show how the effect of de- 
tector efficiency modifies the correlations caused by the 
transverse momentum conservation. As seen, the random 
losses of particles lead to the dilution of the correlation 
that is $4, monotonically goes to zero as the fraction of 
registered particles / reg — > 0. 



VII. PROTON-PROTON COLLISIONS IN 
PYTHIA 

After the discussion of various mechanisms respon- 
sible for azimuthal correlations, let us now consider 
more realistic situation where several mechanisms of 
azimuthal correlations are present at the same time. 
We used the PYTHIA generator [T^] to simulate p-p 
collisions at several collision energies accessible at 
SPS (t/sjv^ = 6.27, 7.62, 8.73, 12.3, 17.3 GeV) 
and RHIC (v^jvjv = 19.6, 62.4, 130, 200 GeV). 
For every energy a set of minimum bias events was 
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FIG. 8: (Color online) $0 as a function of fraction of regis- 
tered particles / rcg . The correlation results from the trans- 
verse momentum conservation. 



collected. We treated as stable the following particles: 
[i-, w°, _vr+, K°, K+, K° L , K° s , A, S+, E-, S°, S", Q" 
and their antiparticles. No acceptance cuts were applied. 

For every energy we computed $0 separately for pos- 
itive, negative and all charged particles. The results are 
shown in Fig. [9] As seen, $0 is negative and weakly de- 
pends on collision energy. To understand why <I>0 is so 
different for negative and for positive particles, we ex- 
cluded protons from all charged and from positive parti- 
cles. The corresponding values of <I>0 are also shown in 
Fig. [9j After excluding protons, the correlations among 
negative particles and among positive are very similar to 
each other. 

What is the mechanism responsible for negative values 
of $0? We first checked that high px particles play no 
role here, as $0 does not significantly change when par- 
ticles with pt > 1.5 GeV are excluded. The correlations 
among charged particles can be caused by the effect res- 
onance decays but the effect is certainly very minor for 
same-sign particles, as there are very a few resonances 
decaying into two positive or two negative particles. 

The transverse momentum conservation, which is dis- 
cussed in Sec. VI, is another possible source of neg- 
ative values of $0. We checked that the PHYTHIA 
events indeed obey the transverse momentum conserva- 
tion. Specifically, we proved vanishing of the total mo- 
mentum in x and in y directions of all particles (charged 
and neutral) from every event. To quantitatively study 
the effect of transverse momentum conservation we pro- 
ceeded as follows. For every collision energy we deter- 
mined the average multiplicity of positive, negative and 
neutral particles. Then, we performed the simple simu- 
lations described in Sec. VI, generating events of a given 
total multiplicity which satisfy the transverse momen- 
tum conservation. Then, a fraction of particles was ran- 
domly eliminated to get the multiplicity of charged, pos- 
itive or negative particles. The values of <I?0 computed 



for such events are shown in Fig. 10 As seen, the val- 
ues of <I>0 for the PYTHIA events agree quite well with 
the results of our toy-model simulations which take into 
account only the effect of transverse momentum conser- 
vation. It is somewhat surprising that the agreement for 
same-sign particles is not much better that that for all 
charged particles. It means that the resonance decays do 
not generate strong correlations in the PYTHIA events. 
We note, however, that the effect of transverse momen- 
tum conservation overshoots the correlations of the like- 
sign particles and it undershoots the correlations of all 
charged particles. The latter results presumably signals 
presence of resonances decaying into pairs of one positive 
and one negative particle. 



VIII. NUCLEUS-NUCLEUS COLLISIONS IN 
HIJING 

We have also performed simulations of nucleus-nucleus 
collisions using the HIJING [TS] event generator. We 
have simulated the collisions of p-p, C-C, Si-Si and Pb- 
Pb at ^snn — 17.3 GeV and $0 has been computed 
separately for positive, negative and all charged particles 
coming from minimum bias events. The results are shown 
in Fig. 11 As seen, <I>0 is almost independent of the mass 
number of colliding nuclei and the values of $0 are very 
close to those found using PYTHIA. It is by no means 
accidental. When a nucleus-nucleus collision is a simple 
superposition of nucleon-nucleon interactions, the value 
of $ is exactly the same for p-p interactions and nucleus- 
nucleus collisions at any centrality. In the HIJING model 
a nucleus-nucleus collision is not exactly a superposition 
of nucleon-nucleon collisions but it is almost so. And 
the treatment of proton-proton interactions is essentially 
the same in PYTHIA and HIJING. For these reasons 
our analysis of PYTHIA events presented in the previous 
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FIG. 9: (Color online) The energy dependence of for 
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simulations of p-p collisions. The pale symbols correspond 
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section fully applies here. 



IX. SUMMARY AND OUTLOOK 

Azimuthal correlations of final state particles from 
high-energy collisions carry valuable information on the 
collision dynamics. It motivates the analysis of experi- 
mental data collected by the NA49 and NA61 Collabora- 
tions which is in progress with some preliminary results 
already published [3]. The integral measure which 
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FIG. 11: The Hijing simulation of the system size dependence 
of for nucleus-nucleus collisions at ^/sjviv = 17.3 GeV with 
no acceptance cuts. 



proved to be very sensitive to various dynamical correla- 
tions, is used in the analysis. To interpret the experimen- 
tal results it should be understood how different sources 
of correlations manifest themselves when measured by 
means of This was the aim of our study. We per- 
formed several simulations to analyze separately the az- 
imuthal correlations caused by the elliptic flow, resonance 
decays, jets and transverse momentum conservation. We 
also discussed how the correlations are diluted due to 
randomly lost particles. Finally we used the PYTHIA 
and HIJING event generators to produce a big sample 
of events which mimic experimental data from p-p and 
nucleus- nucleus collisions at the SPS and RHIC collision 
energies. $0 appeared to be surprisingly independent of 
the collision energy and of the size of colliding systems. 
Applying some kinematic cuts and selection criteria of 
particles, we showed that the azimuthal correlations are 
dominated by rather trivial effect of transverse momen- 
tum conservation which appeared to be almost indepen- 
dent of particle's multiplicity which changes dramatically 
for collision energies and system's sizes under considera- 
tion. 



The experience gathered in the course of this theoreti- 
cal study will be used to better understand experimental 
data. Quantitative analysis of several simple mechanisms 
of azimuthal correlations we discussed will facilitate an 
observation of possible new phenomena like critical fluc- 
tuations at phase boundaries of strongly interacting mat- 
ter or plasma instabilities from the early stage of rela- 
tivistic heavy-ion collisions. 
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to back into two particles. When one particle is emitted 
at the azimuthal angle <pi and < <f>\ < 7r, the second 
particle is emitted at 02 = 0i + 7r. When n < 4>\ < 2tt, 
then 4>2 = 4>\— 7T- Therefore, the two-particle distribution 
of azimuthal angles reads 



Pi 



2tt 
1 



e(7r-0i)tf(0i-0 2 + 7r) (A.2) 



- ©(01 - 7T) <5(0i - 02 ~ 7T) 

Z7T 



Appendix: Toy model of resonance decays 



One observes that 



The inclusive distribution of azimuthal angle is as- 
sumed to be flat that is 



PincW = ^ 6(0)6(277 
Z7T 



which gives = it and 2 = Air 2 /3. Consequently, z 2 — 

7T 2 /3. 

Let us first assume that all produced particles come 
from heavy resonances which are at rest and decay back 



^2 



,0 2 )=iinc(0), (A.3) 



(A.l) and computes 



#1<20201 02^(01, 02) = 77 7T 2 

6 



(A.4) 



We further assume that the particle multiplicity is arbitrary but fixed even number N. Then, the N— particle 
distribution of azimuthal angles is 



-FW(01, 02, ■ ■ • , 4>n) = ^(01, 02) f*2(03, 04) ' ' ' Pi (0JV- 1 , 4>n) ■ 

The variable Z is defined as Z — 0i + $2 + ■ ■ ■ + 4>n — Nir and one computes (Z 2 ) in the following way 



(Z 2 ) = / (20^02 . . . #Ar(<Al + 02 + ■ • • + 0W - N-K) 2 P N {(t>l, 2 , . . . , 0jv) 



A0 2 + N 



2„2 



620i02-P 2 (0i, 2 ) + iV(A^ - 2)0 2 - 2iV>7r + N it 



N 

r 



(A.5) 



(A.6) 



Using the result (A.6 1 and keeping in mind that z 2 — 7r 2 /3, we find the formula ^ from the definition Q. 

Let us now assume that N± particles come from resonances and additional N2 particles are produced independently 
from each other and from resonances. Then, we still have z 2 = 7r 2 /3 and (Z 2 ) is computed as 



(Z 2 ) = N(f> 2 + N 1 



«»201 02^2(01,02) +N 1 (N 1 - 2)0 2 



(A.7) 



+ 2N 1 N 2 (j) 2 + N 2 (N 2 - 1)0 2 - 2A 2 0tt + N 2 tt 2 = < — (2iV — AM 

6 



where N = Ni + N 2 . The result ( |A.7| with / = Ni/N gives the formula 

The model can be easily generalized to the situation when the particles from a correlated pair are not emitted back 
to back but their relative azimuthal angle is fixed and equal A0. Then, the two-particle distribution of azimuthal 
angles is 

P 2 (0i, 2 ) = ^ e(27r - A0 - 00 (5(0! - 02 + A0) + -!- 9(0! - 2tt + A0) (5(0! - 2 + A0 - 2tt) , (A.8) 

Z7T ZTT 

I 



and instead of Eq. (A.4) we have 

J d0ld0 2 0102-P 2 (01,02) = A07T+i(A0) 5 



Substituting the result (A.9) to Eq. (A.6) or Eq. (A.7 1, 



(A.9) 



one finds (Z 2 ) which finally gives the formula (10) or (12 1, respectively. 
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